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the researchers in recent years. This can be seen through the various works reported in the literature. 
This paper reviews the various procedures adopted for the exergy based performance evaluation of LHTS 
units. The influence of operating and design parameters on the exergy stored/retrieved and thus, on the 
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optimization is addressed as a main aspect. The need of exergy analysis for the comparative evaluation of 
LHTS systems with performance enhancement techniques is emphasized. Thermoeconomics methods 

applicable to LHTS systems are also presented in this paper. 
© 2010 Elsevier Ltd. All rights reserved. 
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Nomenclature 

c specific heat (J/kg K) 

Ex exergy (J) 

Ex exergy rate (W) 

f liquid fraction 

h specific enthalpy (J/kg) 

L latent heat (J/kg) 

M mass of PCM (kg) 

m mass flow rate of HTF (kg/s) 
Ns entropy generation number 
NTU number of transfer units 

P pressure (Pa) 

Q instantaneous heat transfer rate (W) 
R ideal gas constant (J/kg K) 
S entropy (J/K) 

T temperature (°C or K) 

Tm melting temperature (°C or K) 
t time (s) 

V volume (m7?) 

Greek symbols 

p density (kg/m?) 

n energy efficiency 

y exergy efficiency 

Subscripts 

char charging 

dis discharging 

final final condition 

gain gained amount 

HTF heat transfer fluid 

in inlet 

init initial condition 

l liquid phase 

loss lost amount 

o environment 

out outlet 

PCM phase change material 

S solid phase 

stored stored amount 

wall heat exchanger wall 


major problem in managing energy from the above-mentioned 
sources is the time gap between the availability and the need. 
Employment of effective energy storage/retrieval device is 
emphasized, if the energy source is intermittent and time 
dependent. For the last three decades, there has been a growing 
interest in latent heat thermal storage (LHTS) technique, which is 
proved as a better engineering option over sensible heat storage. 
The relative merits of LHTS include large energy storage for a given 
volume, uniform energy storage/supply, compactness, etc. 

In LHTS units, phase change materials (PCM) are employed, which 
undergo change of phase (solid to liquid and vice versa) during the 
energy transfer process. In addition to the above stated advantages, 
extensive studies on wide range of PCMs, have made LHTS unit as a 
potential system in various engineering fields. For comprehensive 
details of various PCMs and their applications, readers are referred 
to the recent review articles by Farid et al. [1], Sharma and Sagara 
[2], Zalba et al. [3] and Kenisarin and Mahkamoy [4]. 


However, the employment of LHTS units for energy storage/ 
retrieval can be realized only if the performance of the unit is 
known. The performance assessment of any LHTS unit requires a 
thorough knowledge about the thermal behavior of PCM that is 
employed. Accordingly, a great amount of work both numerically 
and experimentally has been devoted on the study of thermal 
behavior of different configurations of LHTS units applicable to 
various fields. Some of the works are also extended to investigate 
the influence of various geometric, operating and design param- 
eters on the thermal behavior. Investigations by Trp [5] (shell and 
tube LHTS for solar thermal applications), Regin et al. [6] (LHTS 
loaded with cylindrical capsule PCM for solar water heater), 
Zukowski [7] (LHTS with encapsulated PCM for solar space 
heating), Chen et al. [8] (PCM packed in a box-type solar cooker), 
Akgun et al. [9] (PCM packed in a conical shell and tube LHTS) and 
Guo and Zhang [10] (Aluminium foil attached LHTS for solar power 
generation) are some notable works in recent years focusing on the 
above-mentioned points. 

Once the thermal behavior is known, the performance of LHTS 
units can be evaluated through two approaches. They are: 


e By applying energy conservation principle (based on first law of 
thermodynamics). 

e By applying exergy principle (based on second law of thermo- 
dynamics). 


The effective handling of performance evaluation of LHTS 
unit would result in optimized system for the application of 
interest. 

However, in a recent review, Verma et al. [11] have stressed the 
need for exergy analysis for LHTS units. Based on this, it is felt that 
the exergy analysis should be taken as a main aspect in the future 
research works on LHTS units. This paper is motivated from the 
above point and is intended to review the works on various exergy 
based performance evaluation techniques employed for LHTS 
units. This review is intended to provide necessary guidelines for 
the future exergy based studies. 


2. Energy analysis of LHTS systems 


The performance of a thermal system can be evaluated through a 
well known parameter, efficiency or effectiveness. As far as LHTS 
systems are concerned, the efficiency is nothing but a measure of 
how effectively the heat or cold energy is stored or recovered. In 
other words, it indicates how much heat or cold that could not be 
stored or recovered. Hence, the efficiency of LHTS can be determined 
as charging efficiency or discharge efficiency or overall (complete 
cycle) efficiency. The terms ‘charging’ and ‘discharging’ are used, 
respectively, for melting and solidification processes of PCM. This is 
because the heat energy is stored in the PCM during melting and the 
same is retrieved from the PCM during solidification. However, 
when it comes to application like air conditioning systems, the 
desired form of energy is ‘cold’ which implies that the cold energy is 
stored during solidification and it is retrieved during melting. In this 
case ‘charging’ and ‘discharging’ can be applied to solidification and 
melting, respectively. In view of generalization and to avoid 
confusion, however, we adopt ‘charging’ for melting (corresponds 
to heat gained or cold released by PCM) and ‘discharging’ for 
solidification (corresponds to heat released or cold gained by PCM). 

Many investigators carried out calculation of efficiency or a 
similar parameter of LHTS units employed for different applica- 
tions. These investigations have focused on either charging and 
discharging processes separately or the complete working cycle of 
the system depending on the application. A system is always 
expected to produce the maximum possible efficiency, but the 
system with maximum efficiency is possible only at certain 
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conditions. Hence, the influence of various design and operating 
parameters on the efficiency has also been addressed. 

El Qarnia [12] calculated the storage efficiency of LHTS used for 
solar water heater. The efficiency is expressed as a ratio between 
the latent heat stored in the PCM and the total solar radiation. 
Higher mass flow rates of heat transfer fluid (HTF) and more 
number of tubes used in the heat exchanger resulted in maximum 
storage efficiency. In the similar way, Kaygusuz [13] has computed 
the storage efficiency of LHTS unit of a solar heat pump. It is 
reported that stored heat and storage efficiency of LHTS were 
found to be increasing with increase in mass flow rate of HTF. In 
addition to that, the effect of inlet temperature of HTF on the 
storage capacity was studied and the effect is reported as less 
significant. 

Seeniraj et al. [14] investigated the thermal performance of 
shell and tube LHTS unit employed for space based power 
generation. The thermal performance was studied using a quantity 
which is the ratio between the total heat stored and the maximum 
latent heat that can be stored. The results have once again proved 
that higher mass flow rates of HTF would lead to high storage 
performance. Moreover, it is shown that the storage performance is 
higher in case of smaller size unit as compared to that of larger one. 

Although Canbazoglu et al. [15] have not calculated the 
efficiency, the total heat (sensible + latent heat) that can be stored 
in the PCM was calculated to compare the storage capacity of 
various hydrated salts. The LHTS unit taken up for the study is the 
one integrated with solar water heaters. 

In some applications, the performance of LHTS system must be 
known during discharging period rather than during charging 
period. Gumus [16] has designed a LHTS unit for preheating of 
internal combustion engines at cold conditions. Hence, the 
efficiency was calculated during discharging process which is 
expressed as ratio of energy gained by engine components and 
total heat stored in the unit. It is found that the discharge efficiency 
of the unit increases with time and the unit could produce a 
maximum efficiency of 57.5%, which is obviously at the end of the 
discharging process. 

Devahastin and Pitakusuriyarat [17] proposed a shell and tube 
LHTS unit for drying food products like sweet potato. In this case, 
performance during the heat recovery process by the HTF is 
obviously more important. Hence, the authors calculated the rate of 
extractable energy from the unit using the temperature difference 
between the inlet and outlet of the HTF. It is also reported that the 
extractable energy per unit mass flow rate of HTF decreases as the 
velocity of HTF increases. The rate of heat recovered from the unit is 
taken as a performance parameter also by Mettawee and Assassa 
[18]. The heat recovered was calculated from the heat gained by the 
HTF. The effect of mass flow rate on the heat recovery was found to 
similar to that reported in the other works. 

The calculation of efficiency during both charging and 
discharging processes but separately can also be found in some 
works. Kaizawa et al. [19] have used heat storage ratio and heat 
release ratio as performance parameters. The heat storage ratio is 
given as, 


Total energy stored in the unit 
Maximum storage capacity of the unit 


Similarly, the heat release ratio is expressed as, 


Total energy released by the unit 
Maximum storage capacity of the unit 


The heat stored and released are calculated from the heat 
absorbed and released, respectively, by the HTF. The maximum 


heat storage capacity takes into account the sensible and latent 
heat capacity of the PCM. The results show that both the 
performance ratios increase with increase of mass flow rate of 
HTF at all times. 

Mawire and McPherson [20] have defined charging and 
discharge efficiencies separately for a LHTS unit proposed for 
solar cookers as follows: 


_ [ Total energy stored in the unit 
"ehar = | Total energy supplied by the HTF] ’ 
= ae energy absorbed by the | 
dis 


Total energy stored in the unit 


However, the results of these efficiencies are not reported, as 
the objective of the work was limited to development of control 
techniques for control problems. 

The complete cycle or overall efficiency of the system has also 
been used by some researchers. Mettawee and Assassa [21] have 
defined mean daily efficiency for a LHTS integrated with a solar 
collector which is used for hot water production. The mean daily 
efficiency considers total heat absorbed by HTF and total incoming 
radiation as output and input, respectively. Exactly in the similar 
way, Sharma et al. [22] have calculated a performance parameter 
to represent the overall performance of a LHTS of solar cooker. 
However, the effect of operating/design parameters on the 
efficiency is not discussed in these papers. 

From the summary given above, one may think that detailed 
information on system efficiency or capacity either during 
charging and discharging separately or during the complete 
working cycle along with the influence of various operating/ 
design parameters are explored. Subsequently, the explored 
information would enable the designer to develop an optimum 
system. However, the performance of LHTS units is assessed using 
only energy balance in these studies. This procedure, in fact, is 
based on first law of thermodynamics and generally employed to 
identify the ways to improve the quantity of heat stored/recovered. 
At this point, it is important to mention the popular quote by Bejan 
[23] which is 


“the primary purpose of a thermal energy storage system is not, as 
the name implies, to store the energy, rather, to store useful work”. 


From this perspective, it can be stated that first law based 
analysis, i.e. energy analysis becomes inadequate as information 
on usefulness of energy cannot be obtained. In other words, energy 
analysis does not reflect the quality of energy stored/recovered. 
Hence, an ‘optimized LHTS unit’ based on energy analysis is not 
necessarily an optimum system. This inadequacy can be overcome, 
if an analysis based on second law of thermodynamics is carried 
out for LHTS units. This is because the degradation of energy 
stored/recovered is very much appreciated only in the second law 
analysis. Rosen et al. [24] have reported that energy analysis is 
complicated and confusing when it comes to cold thermal storage 
systems those correspond to one employed for building cooling. 
This is because only heat flows are taken into account rather than 
cold flows, in energy analysis. On the other hand, second law based 
analysis can treat both heat and cold which is out of equilibrium 
with the environment, as a valuable commodity. Hence, second 
law analysis is inherently the better option than energy analysis 
especially in case of cold thermal storage systems. Moreover, a 
thermoecomomics analysis of a thermal system is necessarily 
based on second law analysis in order to obtain feasible design and 
operating conditions. This is because only second law analysis 
correctly reflects the economic value of the storage operation [25]. 
Thermoeconomics analysis combines second law analysis and 
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Table 1 
Benefits of second law analysis in LHTS systems. 
S No Benefits 
1 Quality of energy stored/recovered can be measured 
2 Destroyed energy quality as well as lost energy quality 


can be quantified 


3 Causes of internal and external irreversibilities can be identified 
4 Both heat and cold can be treated as valuable commodities 
5 Provides a meaningful base for themoeconomics analysis 


engineering economics. Since the second law analysis takes into 
account the true potential of energy, estimation of cost and 
subsequent optimization of the system can be realized through 
thermoeconomics. Further discussion on thermoeconomics analy- 
sis of LHTS systems is presented in the final section of this paper. 
The benefits of second law analysis are highlighted in Table 1. 


3. Second law analysis-Thermodynamic concepts 


The second law of thermodynamics is very distinct from the 
first law as it introduces new concepts of exergy and entropy. 
Hence, the analysis based on the second law of thermodynamics 
needs to be established with a clear understanding of these 
concepts. This section overviews the concepts of exergy and 
entropy and their roles in the second law based analysis. 


3.1. Exergy 


Unlike first law of thermodynamics, second law of thermody- 
namics quantifies the quality of energy. The quality of energy is 
gauged from the state of system in relation to the surrounding 
conditions. This quality or usefulness of energy is termed as 
“exergy”. Many more equivalent terms like availability, available 
energy, essergy, work capability, utilizable energy, etc., can be 
found in literature. However, Kaygusuz and Ayhan [26] have 
reported that use of the term exergy should be encouraged, as it 
was agreed at the “Fourth International Symposium on Second Law 
Analysis of Thermal System” held at Rome, Italy in 1987. Moreover, 
the term exergy is more commonly employed than others. Hence, 
in this paper, we are confined to the term exergy. Since second law 
of thermodynamics introduces “exergy”, the analysis based on 
second law is widely called exergy analysis. 

Exergy is defined as maximum quantity of work that can be 
produced by a system as it comes to equilibrium with surrounding. 
In other words, exergy is nothing but a measure of potential of the 
system to cause change, due to its non-equilibrium condition with 
surrounding. A system can carry exergy only ifitis notin equilibrium 
with surrounding and the exergy becomes more and more as the 
system deviates more and more from surrounding. Once the 
equilibrium is attained then exergy becomes zero. This means 
exergy cannot be conserved, but can be consumed or destroyed. This 
statement highlights the reason why energy is getting degraded. 

It is clear from the abovementioned points that the degradation 
of energy is automatically accounted for in the analysis, if it is 
second law (exergy) based analysis. As mentioned earlier, the LHTS 
system is for storing useful work. Also the recovery of useful work 
is equally important. Hence, the calculation of exergy transfer is 
more appropriate than the calculation of energy transfer during 
charging and discharging processes. Since a LHTS unit is concerned 
with only heat transfer, the exergy depends on the temperature at 
which it is available in relation to the temperature of surrounding. 


3.2. Entropy 


The exergy content of a system or a matter can be obtained 
using the thermodynamic relations known as entropy relations. As 


the name implies, these relations are basically expressed for a 
property known as entropy. The association between entropy and 
second law of thermodynamics is well known as the latter is often 
called law of entropy. 

Moreover, as we have seen before, unlike energy, exergy cannot 
be conserved, rather is destroyed. However, exergy can be 
conserved, if the process takes place in a reversible manner. It is 
well known that all real processes are irreversible and hence, the 
irreversibilities associated with the process can be stated as 
responsible for the exergy destruction. The exergy analysis is not 
only concerned with exergy, but also with exergy destruction. This 
means that it is also necessary to know that to what extent exergy 
can be destroyed in order to quantify the true potential of the 
system/process. This is where the concept of entropy becomes 
more useful. 

Entropy here is used as a measure of irreversibility. This is 
explained as follows. Any process which is irreversible, is 
accompanied by energy degradation and when the equilibrium 
is attained, the exergy content reaches zero value. At the same 
time, the irreversible process is also accompanied by increase in 
entropy which reaches a maximum when the system comes in 
equilibrium with surrounding. Thus, all irreversible processes 
proceed in the direction of increase in entropy and decrease in 
exergy. This highlights the relation between exergy and entropy. In 
actual system, entropy is always generated and exergy is 
destroyed. Hence, the destroyed exergy must be proportional to 
generated entropy. As standard entropy relations are much easier 
to deal with, the generated entropy is included in the exergy 
analysis to explore the destroyed exergy quantity. The generated 
entropy is expressed in terms of entropy generation number which 
will be discussed further in the subsequent sections. 


3.3. Second law efficiency 


First law based efficiency can be stated as the ratio of energy 
output and energy input. In practical thermal systems the energy 
output is less than energy input because of energy loss. The first 
law efficiency indicates the amount of loss that occurs during the 
process in the system. Hence, efficiency can be improved only by 
reducing losses. However, the losses do not reflect the degradation 
of energy. To LHTS units, this is of more important, which can be 
explained as follows: 

The efficiency of a well insulated LHTS unit during the cycle 
(comprises of melting and solidification processes) can be defined 
as the ratio between total energy recovered during solidification 
and total energy stored during melting. Since the latent heats of 
both melting and solidification are the same, the efficiency of the 
LHTS unit becomes 100%. Due to the adiabatic condition, no energy 
loss can occur. Nevertheless, energy would be destroyed due to 
internal irrerversibilities. Even if the realistic condition of non- 
adiabatic is included in the analysis, the energy efficiency cannot 
reflect the destroyed quantity as it measures only energy loss due 
to infiltration. On the other hand, the second law or exergy 
efficiency which takes into account the exergy, can identify and 
quantify both the destroyed quantity as well as lost quantity. Due 
to this reason, the exergy efficiency is found to be less in LHTS 
systems when compared to energy efficiency [24,27-34]. This can 
be seen from Table 2 in which the investigations presenting the 
comparison between energy and exergy efficiency are summa- 
rized. 

The second law or exergy efficiency of a thermal system is thus, 
defined as the ratio of exergy output and exergy input. For a given 
quantity of exergy, the output exergy is less due to exergy 
destruction. As we know exergy destruction is due to irrerversi- 
bilities. Thus, exergy efficiency is a measure of irrerversibilities and 
is defined as, 


2584 


Table 2 
List of works presenting the comparison between energy and exergy efficiencies. 
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Ref. LHTS module Application Operation mode Exergy efficiency based on Exergy efficiency in comparison 
with energy efficiency? 

[24] Rectangular tank Cold thermal storage Cycle Exergy contents Less by 80% 

[27] Not available General Cycle Exergy contents Less by 60% 

[28] Cylindrical tubes Solar space heating Cycle Exergy contents Less by 37% 

[29] Cylindrical tubes Solar space heating Cycle Exergy contents Less by 35% 

[30] Cylindrical tank Solar green house Charging Exergy contents Less by 25% 

[31] Rectangular tank Solar collector Charging Entropy generation Less by 95% 

[32] Spherical capsules General Charging Entropy generation Less by 83% 

[33] Shell and tube General Discharging Exergy contents Less by 47% 

[34] Spherical capsules Cold thermal storage Charging and discharging Entropy generation Less by 7% (charging); Less 


* Approximate representative values. 


Exergy destroyed 
Exergy input l 

As already mentioned, the more exergy is destroyed the more is 

the generated entropy. Hence, the ratio of destroyed exergy and 

exergy input can be replaced by a factor called entropy generation 

number (N,). This implies that, 


w=1-N, (1) 


Eq. (1) shows that the system performance can be improved by 
minimizing the entropy generation number. The entropy minimi- 
zation is taken as an important criterion in the thermal system 
optimization [35]. 

From the above discussions, it is clear that exergy and entropy 
are closely associated with each other and the analysis based on 
second law of thermodynamics revolves around these two. For 
more comprehensive discussion on exergy and entropy, it is 
suggested to refer to the paper by Dincer and Cengel [36]. 


4. Exergy analysis of LHTS systems 


In general, exergy analysis is based on the thermodynamic 
concepts associated with second law of thermodynamics along 
with conservations principles of mass and energy. It is widely 
proved that exergy analysis is very effective tool or compliment to 
energy analysis for performance assessment and optimization of 
thermal systems. This section presents the exergy methodologies, 
i.e. methods of evaluation of exergy based performance parameters 
applicable to LHTS systems. 


4.1. Exergy efficiency 


Rosen [37] has pointed out that though energy based 
efficiencies for thermal energy storage systems are reasonable 
and widely applied, more meaningful efficiencies are defined 
based on exergy. 

The exergy efficiency is evaluated by establishing the exergy 
analysis and subsequently, the system can be optimized in sucha 
way that the exergy efficiency of the system is maximum possible. 
Similar to energy efficiency, the exergy efficiency is expressed for 
charging period and discharging period separately or for the 
complete cycle. 

During charging mode, the HTF transfers the exergy to PCM and 
part of the exergy is stored in the PCM. Hence, the exergy stored is 
considered as desired output in the definition of exergy efficiency. 
Accordingly, the exergy efficiency is expressed as given by 
Watanabe and Kanzawa [38], 


by 14% (discharging) 


Exergy stored in the PCM 


Wenar = Exergy supplied by the HTF 


Since in LHTS systems the heat transfer is time dependent, it is 
also important to evaluate the exergy efficiency at different times 
during melting. This prompts to define the exergy efficiency in 
terms of exergy rate. So, 


Rate of exergy stored in the PCM 
Rate of exergy supplied by the HTF 


Wehar = 


However, Rosen and Dincer [39] have reported that the power 
input to handle the HTF (pump work) must be accounted for in the 
exergy efficiency evaluation. The authors have demonstrated that 
the difference between the exergy efficiency which neglects the 
pump power and that of considering pump work is significant. 
Moreover, this difference is more pronounced in case of exergy 
efficiency than in case of energy efficiency. An expression for 
exergy efficiency considering the pump work can be found in [30] 
and the same is given below: 


Rate of exergy stored in the PCM 


Rate of exergy supplied by the HTF 
+Power input to compressor or pump 


Wehar a 


In the above expressions, the exergy supplied by HTF is nothing 
but the change in flow exergy of HTF during heat transfer with 
PCM. The change in flow exergy can be calculated by using the 
difference between the inlet and outlet temperatures of HTF along 
with the environment temperature. Hence, the exergy efficiency 
indicates the amount of exergy stored from the supplied exergy 
during the heat transfer. The maximum exergy that can be 
supplied by the HTF, however, is not considered here and hence, 
the maximum possible exergy storage capacity of LHTS system 
cannot be obtained. This deficit can be overcome by evaluating the 
exergy efficiency as defined by Gong and Majumdar [40] and the 
same is given as, 


Rate of exergy stored in the PCM 


Venar = Exergy rate possesed by the HTF before contact with PCM 


Similar expression is also reported by Demirel and Ozturk [41]. 
It is clear from the expression that this form of exergy efficiency 
takes into account the temperature at which the exergy is available 
for storing as it uses the difference between the inlet temperature 
of HTF and environment temperature rather than the difference 
between the inlet and outlet temperatures of HTF. Thus, it 
obviously provides the information on how much maximum 
exergy that can be stored from the maximum available potential 
work. 
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Similar expressions can be obtained also for discharging 
(solidification) mode. In the discharging mode, the exergy output 
is the exergy gained by the HTF which is the change of flow exergy 
of HTF. On the other hand, the exergy input is exergy available with 
the PCM. Hence, the following expression can be given for the 
exergy efficiency of the system during the discharging process [38]. 


Exergy gained by the HTF 
Initial exergy available with the PCM 


Wais _ 


This expression provides the total exergy extracted by the HTF 
from the maximum available exergy with the PCM. 

From the perspective of time dependent operation, the exergy 
efficiency can be defined as [27,40], 


Rate of exergy gained by the HTF 
Rate of exergy released by the PCM 


Wais _ 


One can identify the variation of exergy efficiency with respect 
to time during discharging process from the above expression. 

The operation of thermal energy storage systems is inherently a 
cycle comprises of energy storage process followed by energy 
removal process. The analysis may lead to serious errors, if this 
cyclic nature is not accounted for, i.e. analyzing either charging or 
discharging alone [42]. Although the author has highlighted this 
point with reference to sensible heat storage systems, it is felt that 
the same is applicable to LHTS systems too. This is demonstrated 
by Bellecci and Conti [43] as they could minimize the errors arising 
from the latter approach by switching over to the former. The study 
focused on LHTS integrated solar assisted heat engine, nevertheless 
the overall exergy efficiency calculation is recommended for all 
applications involving LHTS units. This can also be justified from 
the results of Gong and Mujumdar [40]. The authors have proved 
that analysis considering charging alone underpredicts the exergy 
efficiency significantly. Hence, the calculation of overall cycle 
exergy efficiency becomes necessary. 

The evaluation of exergy efficiency can readily be done from the 
exergy efficiencies of charging and discharging processes. The 
required expression is as follows [27], 


Woverall = Wehar i Wais (2) 


Alternatively, the overall efficiency is given as [28,29], 


Exergy extracted from the PCM by the HTF during discharging 


Woveratt = Exergy input to the PCM during charging 


Itis evident from the above discussions that the several forms of 
exergy efficiencies are developed and used by the investigators. 
The various forms of exergy efficiencies are summarized in Table 3. 
Nevertheless, the evaluation of overall exergy efficiency is proved 


to be better approach than limiting to charging or discharging 
modes alone. 


4.2. Exergy evaluation 


The previous section illustrated that the evaluation of exergy 
efficiency for LHTS systems requires evaluation of exergy 
associated with PCM and HTF. This section addresses the issues 
related to calculation of exergies of LHTS systems. Demirel and 
Ozturk [41] have employed the following expression for the rate of 
exergy supplied by HTF during charging period. 


: f Teas 
EXin put = MureCure | (THrF in — Torr our) — To In (7m) | (3a) 
HTF out 


Table 3 
Various forms of exergy efficiencies. 
Efficiency Expression Description Ref. 
Charging (Wenar) 1. eet Presents total exergy 38] 
i stored out of supplied 
2 Pamei Presents time-wise 26] 
ai variation of exergy 
f efficiency 
Bs eS Takes into account 30] 
XyrF +Pum p = work n 
t the pumping power 
4. uL Presents maximum 40,41] 
XHTF init A 
possible exergy stored 
Discharging (Wais) 1. re Presents total exergy _ 38] 
. recovered out of supplied 
2 pue Presents maximum 27,40] 
XPCM . 
possible exergy recovered 
Overall (Woveralt) il oe Presents total exergy 27] 
recovered out of supplied 
Ds Wehar Wais mat 28,29] 
Charging/ 1.1—N, Presents the quantity 52] 
discharging/ of exergy destroyed 


overall 


Eq. (3a) expresses exergy input as the change in flow exergy of 
HTF as a result of heat transfer. However, as stated earlier, the 
exergy input can be the exergy content of HTF. The exergy content 
of the HTF is the minimum useful energy obtainable as it reaches to 
the state of environment. Therefore, the rate of exergy input can be 
expressed as given in Eq. (3b). 


F ; Tes 
EXin put = Mate Cute | (THTF in — To) — To In ( me) (3b) 
oO 


As it is known, the total exergy input to the PCM during any 
time interval can be readily obtained by the time integration of the 
relevant expressions. 

During the charging mode of operation, the exergy stored at any 
time instant is computed from the instantaneous heat transfer rate. 
The instantaneous heat transfer rate (Q) can be obtained by 
establishing the energy balance i.e. 


Heat gained by PCM = Heat trasnferred by HTF 
Thus, 
Q = Mute Cure (Tye in — THTF our) (4) 


Now, the rate of exergy stored in PCM is given as, 


EXstored = aft = ( To J| (5a) 


Tpcm 


The temperature of the PCM (Tpcm) is taken as its melting 
temperature by Gong and Mujumdar [40], as the PCM is assumed 
to be at its melting temperature right through the melting process. 
This means the sensible heat of PCM is neglected. However, Farid 
and Kanzawa [44] have stated that in practical LHTS systems, the 
sensible heat contributes considerably to the total heat. This is 
because the PCM is subjected to sensible heating prior to melting 
and subcooling after solidification. This indicates that the PCM 
temperature may be different from its melting point before and 
after the charging/discharging process. Fig. 1 illustrates the 
temperature profile of PCM during phase change with sensible 
heating/cooling. As shown in Fig. 1, the beginning of charging 
process is at its subcooled solid temperature and the charging is 
terminated when the PCM reaches a superheated temperature. In 
this case, the temperature of PCM is computed as the average of 
initial and final temperatures of PCM [41]. However, this is simply 
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T pcm, superheated 


Charging 


ee 


ee 


Discharging 


aS Tpom,subcooled 


Fig. 1. Temperature profile of PCM during phase change with sensible heating/ 
cooling. 


an approximation as its fails to take into account the instantaneous 
temperatures of PCM. Moreover, at any time, the temperature of 
PCM is different at different locations because of the continuous 
movement of solid/liquid interface as the time progresses. Due to 
this fact, it becomes necessary to consider the temperature 
variation in the PCM. Hence, it is suggested to evaluate the Tpcy as 
the average of several temperature values measured/calculated at 
any instant at different locations. 

The energy balance used to obtain Eq. (4) assumes that there is 
no heat loss from the system during the heat transfer, i.e. the 
system is perfectly insulated. Following this, the exergy balance 
can be arrived as, 


Exergy input = Exergy output + Exergy destroyed 


The exergy efficiency based on the above exergy balance can 
reflect only the destroyed exergy due to internal irreversibilities. 
According to Erek and Dincer [33], the heat loss/heat gain is the 
primary parameter in thermal behavior of the system and both the 
energy and exergy efficiencies depend on the heat loss/heat gain. 
Moreover, the authors have reported that better agreement could 
be achieved between numerical and experimental results of energy 
and exergy analysis when the heat gain was taken into account. 
From this perspective, the energy balance for charging process may 
be established as follows in order to include the effect of heat loss/ 
heat gain. 


Energy transferred by HTF = Energy gained by PCM 


+ Energy lost to the surrounding 


Subsequently, the exergy balance is written as, 


Exergy associated with HTF = Exergy stored + Exergy destroyed 


+ Exergy lost to the surrounding 


Hence, the rate of exergy stored is computed from the 
expression given below: 


buns GE) 


where Qioss is the rate of heat loss from the LHTS system due to 

temperature difference between the system and surrounding. 
The total exergy stored in PCM during charging can be readily 

obtained by performing time integration of Eq. (5b). However, it 


requires the computation of heat loss to the surrounding. An 
alternate form, simpler than the time integrated form of Eq. (5b) 
can be found in Ref. [34] as, 


T 
EXstorea = ML i = (7) 
m 


T 
+ Mcpcms (Tn — Tpcm.init) — To In G m | 
PCM. init 


m 


F 2 
+ Mcpcm Toc, Tm) — To in ( refine! (6) 


The first term on the right hand side of Eq. (6) is the exergy 
stored during melting and the second and third terms indicate 
exergy stored during sensible heating before and after the melting, 
respectively. Eq. (6) is not only simpler but also more direct and 
appropriate as it also includes exergy stored during sensible heat 
transfer before and after phase change process. 

Similarly, the more appropriate form of exergy balance for 
discharging process would be, 


Exergy of PCM + Exergy gain from surrounding 
= Exergy gained by HTF + Exergy destroyed 


The exergy balance for a typical LHTS system is illustrated in 
Fig. 2. Based on exergy balance, the corresponding exergy 
expressions for discharging process can be obtained in the similar 
way corresponds to the charging process. The relevant expressions 
proposed by various investigators are listed in Table 4. 


4.3. Entropy generation 


The exergy efficiency computed in terms of exergy output and 
exergy input, for instance, exergy stored in the PCM and exergy 
supplied by the HTF provides the information on how much exergy 
is destroyed. Since this class of expressions does not stress the 
calculation of exergy destruction, one may not be able to identify 
the reasons for exergy destruction. In other words, since the exergy 
destruction is due to the internal irreversibilities, no information 
on when and where the irreversibilities occur, can be obtained. In 
the process of optimization, however, it is the irreversibilities that 


(a) Ex,,,(to surroundings) 
Ex. t + 
EXurF,in m To ExXstored 
= eid + 
Exyrrout 4m 
i 
(b) Ex,,,(from surroundings) 
E EXinput + 
XHTF,out <a Expo 
a quae D 
EXyrF,in m 
oo 
Exgestroyed3, £ 


Fig. 2. Exergy balance for a typical LHTS system (a) during charging and (b) during 
discharging. 
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Table 4 
Expressions of exergy input and output for discharging process. 


Expression 


To 


Exergy input 1. Myre Cue (Ture out — THTF in) [1 = (22) | 


2. Myre Cute (THTF out — Trenein)|1 = ( z )| + Qeain [1 as ( = )] 


Tpcm 


3. miyrrCHrr(Tpcm — Tuten) [1 = (2) 


4.ML|1 = (72)| + Mcpcms [ (Trem ine —Tm) —To In (ae) | + Mepcm [Tm Tpcm, finat) — To In ( = )| 


m 


Exergy output 1. Myre Cure [ (Ta out —Turrin) — To In (zw) ] 


HTF, in 


Description Ref. 
Presents rate of exergy [27] 
released by PCM 

Takes into account the [33] 
heat gain from surroundings 
Presents rate of [40] 


maximum possible 

exergy released by PCM 

Presents total exergy [34] 
released by PCM 
Presents rate of 

exergy recovered by HTF 


Tpcm. final 


[27,40] 


are to be identified so that an attempt to reduce the same can be 
made. 

The difference between reversible and irreversible process is 
described quantitatively by the amount of entropy generation. 
Thus, the entropy generation is the direct measure of internal 
irreversibility or exergy destroyed. The optimum system can now 
be developed in which the entropy generation is minimal. In fact 
Bejan [35] has stated that optimization of a system through 
entropy generation minimization is distinct from that of through 
exergy analysis. This is basically because the exergy analysis deals 
with only thermodynamics where as to calculate and minimize 
entropy generation one needs to use heat transfer and fluid 
mechanics principles along with thermodynamics. 

For LHTS systems, during charging or discharging process, the 
net entropy generation (Sgen) is written as [32,45-47], 


Sgen = S flow + ASpem + ASute (7a) 


where Sjow is the entropy transfer due to the heat transfer between 
PCM and HTF. The second and third terms on the right hand side of 
Eq. (7a) correspond to the entropy variations in the PCM and HTF, 
respectively. 

In the works of Refs. [32,45-47], the PCM was encapsulated in 
plastic shells and HTF was made flowing through the tank filled 
with PCM capsules. Hence, the entropy variation of plastic capsules 
could be neglected. However, other configurations of heat 
exchangers are also investigated for LHTS units employed for 
various applications like solar heating systems, A/C systems, waste 
heat recovery systems, space based solar power systems, etc., For 
the recent review on various kinds of heat exchangers proposed by 
researchers, the readers are referred to [48]. Among these shell and 
tube/concentric double pipe heat exchangers are proved as high 
efficient for minimum volume [49]. In shell and tube configuration, 
the wall of the heat exchanger plays a crucial role in the heat 
transfer as Jegadheeswaran and Pohekar [50] have proved that 
increasing the wall thickness leads to lesser heat transfer rate. 
Therefore, the entropy variation of wall needs to be included in the 
computation of net entropy generation. Accordingly, Erek and 
Dincer [51] have suggested an expression for entropy variation, 
which further includes the entropy generation due to the 
interaction of system with the environment (Shear toss/gain). Since 
the authors investigated solidification of PCM, entropy generation 
due to heat gain by the system from the surroundings is included. 
The general form of the expression applicable to both charging and 
discharging processes can be, 


OF Qeai 
Sgen = S pow + ASpcm + ASute + ASwanr + (= T Sas) (7b) 
oO 


In Eq. (7b), the positive and negative signs stand for Qjoss and 
Qeain, respectively. 


The thermal storage is necessary under the following two 
circumstances: 


e Heat source is available excessive than the demand, so the 
excessive energy would be stored. 

e Heat source is available but there is no demand, thus, entire 
energy is to be stored. 


In the former case, after the storage, the HTF is used in the 
application. Therefore, there will be no entropy generation due to 
the difference between HTF outlet temperature and surrounding 
temperature (Fig. 3a). On the other hand, the latter case requires 
discharge of HTF into the atmosphere (Fig. 3b). Due to this, there 
exits an entropy generation and the same has to be included in the 
net entropy generation. Thus, 


or i 
Sgen = S pow 4 ASpcm + ASurr + ASwan + (= T sa) + ASenvn 
(6) 


(7c) 


The last term on the right hand side of Eq. (7c) denotes the 
entropy generation due to the discharge of HTF into the 
atmosphere. In case of discharging process, HTF will always be 


(a) Qioss 
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Fig. 3. Entropy generation in a typical LHTS system (a) charging—HTF is used in the 
application after storage, (b) charging—HTF is dumped into atmosphere and (c) 
discharging. 
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used in the application after the energy recovery from the LHTS 
unit (Fig. 3c). Hence, ASenvn is always eliminated for discharging. 

The terms on the right hand side of Eq. (7c) are given by the 
following Eqs. 


; £ T 
S pow = rivecine | In (faa) dt (8a) 
0 HTF in 
Eq. (8a), however, ignores the pressure drop irreversibilities due 
to the flow of HTF. To include the same, Kousksou et al. [32] have 
proposed the following expression. 


i T 
S fow = recone |, In (4) dt 


Tut in 


t R 
+ riurR f In (7m) dt (8b) 
0 Purr out 


The second term on the right hand side of Eq. (8b) presents the 
pressure drop irreversibilities. However, the HTF is assumed to be 
an ideal gas, which may not be valid in case of liquids. MacPess and 
Dincer [34] have stated that the entropy generation due to the 
pressure drop irreversibilities is the entropy generation due to 
viscous dissipation. Accordingly, Eq. (8b) can be modified as, 


£ F 
S flow = riecone |, In (4) dt 


Tut in 


_ Myre At(Prre in — PHTF out) 
T 
PHrFTHTF 


(8c) 


where Tyr is the average of initial and final temperatures of HTF 
and the second term on the right hand side of Eq. (8c) is the entropy 
generation due to the pressure drop irreversibilities. 

Now, we address the computation of entropy variation in the 
PCM. El-Dessouky and Al-Juwayhel [52] have assumed that during 
the phase change process, the specific Gibbs free energies for solid 
and liquid phases are same at the melting point. This implies that, 


hpcm,s — TmSpcm.s = pcm — TmSpcm.t (9a) 
hence, 


_ hpcms — Mpem ML 
Tm Tm 


ASpcm (9b) 

Eq. (9b) is valid only if the temperature of PCM is equal to its 
melting point throughout the phase change. For a process in which 
sensible heat transfer is involved, the expression should take the 
following form [47]. 


t 
ASpcm = Ppcm LV Pem | ft) dt 


t 
Tpcm(t) ) 
+ CecmV. 1 In (ee dt 9c 
Ppcm€Pcm V PCM k Trm(t = 0) (9c) 


The first term on the right hand side of Eq. (9c) denotes the 
entropy variation in the PCM due to phase change and the second 


Table 5 
Studies on evaluation of entropy generation. 


term denotes the entropy variation in the PCM due to sensible 
heating/cooling. 

Similarly, the entropy variation in HTF and wall can be written, 
respectively, 


: Ture (£) 
A Sae pics, [ny 10 
HTF = PHTFCHTF V HTF A Tar(t = 0) (10) 
t 
_ Twan(t) 
ASwatt = PwallCwall V wall o In Ce: _ 0) dt (11) 


The literature review reveals that none of the works considered 
all the terms appearing in Eq. (7c) together, for the calculation of 
net entropy generation. This can be seen from Table 5, in which 
various investigations on evaluation of entropy generation are 
listed. However, it is recommended that all the terms appearing in 
the Eq. (7c) need to be considered as each term contributes 
significantly to the net entropy generation. 

Once the entropy generation is known, the destroyed exergy 
(Exdestroyea) Can be obtained from the well known Gousy-Stodala 
theorem [35], which is in the following mathematical form, 


EXdestroyed = ToSgen (12) 
Then the exergy efficiency becomes, 


= on oe] 
4 1 E (13) 

As stated earlier, the ratio of destroyed exergy (ToSgen) and 
exergy input can be replaced by entropy generation number and 
the expression for exergy efficiency takes the form given in Eq. (1). 
It is clear from Eq. (1) that the system achieves the maximum 
efficiency if entropy generation number becomes zero. In other 
words, the system operates in complete reversible manner and 
thus, no exergy is destroyed. This ideal condition is possible only if 
the overall number of transfer units becomes infinite and pressure 
drop is zero [52]. This highlights the importance of entropy 
generation number in optimization of the system. 


5. Exergy based optimization of LHTS systems 


The selection of appropriate values of operating and design 
parameters is a challenging task as different parameters have 
different influence on exergy performance of LHTS systems. This 
section reviews the outcomes of research works reported related to 
exergy based optimization of operating and design conditions of 
LHTS systems. 

In general, one may expect the performance of LHTS systems is 
influenced by the following parameters. 


e HTF 
o Inlet temperature 
o Mass flow rate/velocity 
o Material 


Ref. LHTS module Operation mode Entropy generation terms considered Pressure drop irreversibility 
[32] Spherical capsules Charging S fow:ASpcm, ASHrF Included 
[33] Spherical capsules Charging and discharging S fow:ASpcm, ASue, ASwatt Included 
[45] Spherical capsules Charging S fowASpcm, AStrre Neglected 
[46] Spherical capsules Discharging S fowASpcm, AStrrr Neglected 
[47] Spherical capsules Charging S pow ASpcm, ASpre Neglected 
[51] Shell and tube Discharging S pow ASpem, ASutF, ASwatts Sheat gain Neglected 
[52] Not available Cycle S pow ASpem AStrrASenvn Included 
[55] Shell and tube Charging S pow ASpcm ASyre, ‘ASenun Included 
[58] Shell and tube Cycle S pow ASpcm, ASytr Included 
[77] Not available Cycle S pow ASpem ASprre.ASenvn Included 
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e PCM 
o Melting temperature 
o Initial temperature 
o Material 
e System 
o Dimensions 
o Material 


5.1. Influence of HTF 


The inlet temperature of HTF is very crucial in initiating the 
charging/discharging process. If the temperature of PCM (initial/ 
melting temperature) is fixed, then there exists a question: what 
should be the inlet temperature of HTF? The inlet temperature of 
HTF depends on the heat source during the heat storage and on 
surroundings during heat recovery. However, it should be 
controlled for maximizing the exergy performance. 

During charging, higher HTF inlet temperature leads to higher 
temperature difference between HTF and PCM. This would result in 
higher entropy generation. This is attributed to the fact that the 
entropy generation is directly proportional to finite temperature 
difference. This is proved by El-Dessouky and Al-Juwayhel [52] in 
shell and tube arrangement and more recently by Kousksou et al. 
[32] in PCM capsules system. Hence, one may conclude that for 
minimum entropy generation and thus, for maximum exergy 
performance, the HTF inlet temperature should be as low as 
possible, i.e. close to the initial temperature of the PCM. However, 
the smaller temperature difference would effectively affect the 
heat transfer rate and thus, the charging rate. In LHTS systems 
employed for applications like waste heat recovery systems, the 
charging rate should be high, as the source availability is 
intermittent [53]. Hence, the selection of HTF inlet temperature 
may be a compromise between the required charging rate and 
exergy performance. On the other hand, if the LHTS system is 
employed for solar thermal applications, higher charging rate may 
not be important. This allows lowest possible HTF inlet tempera- 
ture and thus, highest exergy efficiency can be realized. At the 
beginning of discharging process, the PCM in which the energy is 
already stored would be at a high temperature. From the similar 
prospective discussed already, the finite temperature difference 
between PCM and HTF should be reduced in order to decrease the 
entropy generation and exergy destruction. The results of El- 
Dessouky and Al-Juwayhel [52] have proved that entropy 
generation number could be reduced by increasing the HTF (air/ 
water) inlet temperature during discharging, although the effect 
was marginal when water was used as HTF. This may lead to a 
conclusion that the effect of HTF inlet temperature on entropy 
generation is more pronounced when a gas is the HTF than when 
liquid as HTF. However, this general conclusion is not always valid 
as Erek and Dincer [51] have shown significant influence of HTF 
inlet temperature on entropy generation number with ethylene 
glycol as HTF. Unlike El-Dessouky and Al-Juwayhel [52], Venka- 
taramayya and Ramesh [27] and Erek and Dincer [33] have 
calculated exergy efficiency based on exergy contents. In both the 
works carried out with shell and tube module, the authors have 
found significant increase in exergy efficiency due to increase in 
HTF inlet temperature during discharging. Like in charging process, 
operating HTF at a temperature close to PCM temperature is also 
tricky in discharging process, if higher discharging rate is required. 

The exergy performance of LHTS system is also found to be 
dependent on mass flow rate/velocity of HTF. To study the effect of 
mass flow rate/velocity of a fluid, the well known dimensionless 
number called Reynolds (Re) number is generally used. Higher Re 
number indicates higher mass flow rate/velocity and vice versa. 
Recently, Kousksou et al. [47] have reported that the velocity of 


HTF has almost negligible effect on the entropy generation 
number. Since the numerical model neglected the entropy 
generation due to pressure drop in the system, the effects of 
change in velocity (responsible for pressure drop) on entropy 
generation number could not be explored. The increase in Re 
number (higher mass flow rate/velocity) obviously results in 
higher pressure drop and thus, higher entropy generation number. 
However, the results of Kousksou et al. [32] and Erek and Dincer 
[51] and have showed the opposite effect. Erek and Dincer [33] 
have also presented the same results but for exergy efficiency. The 
reason for this can be given as follows. At higher Re numbers, the 
residing time of HTF in the system becomes less which results in 
smaller change in the temperature when HTF comes out of the 
system. Hence, the temperature difference between HTF inlet and 
outlet decreases as Re number increases. Referring to Eqs. (8a) and 
(8b), it can be stated that the smaller temperature difference would 
reduce the entropy transfer due to heat transfer between HTF and 
PCM. Therefore, the net entropy generation can be reduced by 
increasing Re number, although there is an increase in pressure 
drop irreversibility. 

It should be mentioned that none of the above stated works 
included the pumping power in the exergy analysis. It is evident 
that high Re number requires high pumping power which leads to 
increase in mechanical irreversibility. Hence, it is also required to 
take into account the pumping power requirement corresponding 
to Re number, in order to explore the exact influence of Re number 
on entropy generation number. Nevertheless, the influence of Re 
number on exergy performance is reported as less as compared to 
that of HTF inlet temperature. 


5.2. Influence of PCM 


The selection of appropriate PCM for any application is also 
crucial for the proper operation of LHTS systems. In doing so, the 
melting temperature of PCM is taken as a main criterion although 
other factors like latent heat, thermal conductivity, stability, 
availability cost, etc., are taken into consideration. A list of widely 
studied/potential PCMs of different temperature ranges proposed 
for different applications can be found in Ref. [54]. For any 
application, the melting point can be within the operating 
temperature range of the application. Now the question is what 
should be optimum melting point within the range? This is where 
the exergy analysis is of greater help. Kousksou et al. [47] have 
observed that the irreversibility in the LHTS system is strongly 
influenced by the melting temperature of PCM. 

In general, the optimum melting temperature can be obtained 
through a simple procedure as explained by Bejan [35]. The 
operating cycle of LHTS unit consists of heat storage and removal 
processes, can be assumed to follow Carnot cycle. Accordingly, the 
heat is supplied at a constant temperature (melting point) and the 
heat is rejected into a low temperature reservoir (surroundings). 
This assumption leads to formulation of a reversible cycle executed 
by the working fluid working between melting temperature and 
atmospheric temperature. Now, the rate of exergy extracted from 
the PCM is given as, 


; : T 
Ex = MypeCute (Turin — Tm) (1 — e™™Y) (1 = T) (14) 
m 


Maximizing this exergy with respect to melting temperature 
leads to, 


Tm(o ptimum) = \TateinTo 


Hence, the optimum melting temperature is the geometric 
mean of HTF inlet temperature and temperature of surroundings. 
Same condition for optimum melting temperature can be obtained 


(15a) 
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by minimizing the entropy generation [55] or by maximizing 
exergy efficiency [40,56]. Although Eq. (15a) has been arrived 
based on the assumption that PCM is at its melting temperature 
throughout the cycle, according to Bejan [35] it is still a good 
approximation even when sensible heat transfer is involved. 

Recently, Aghbalou et al. [57] have successfully used Eq. (15a) in 
their study carried out on LHTS unit consists of PCM filled 
rectangular slabs. The results showed minimum irreversibility 
when PCM of optimum melting temperature was used. However, 
Charach [58] has pointed out that Eq. (15a) is valid only for 
charging process. As it is known, during the discharging process, 
the entropy generation would be less if the melting temperature 
approaches the HTF inlet temperature. Due to this fact, the 
optimization of melting temperature for the complete cycle 
(charging and discharging) would result in, 
T (0 ptimum) = È (Tyne in + To) (15b) 

Same expression is also presented by Gong and Mujumdar [40] 
which was obtained by maximizing the overall exergy efficiency. 
Charach [58] has also added that the optimum melting tempera- 
ture is close to arithmetic mean of HTF inlet temperature and 
temperature of surroundings only for lower HTF inlet tempera- 
tures. For higher HTF inlet temperatures, the same was found to be 
between the arithmetic and geometric means of HTF inlet 
temperature and temperature of surroundings. On contrary, Conti 
et al. [59] found two distinct values for optimum melting 
temperature on the two sides of (1/2)(Turrin + To) and the lower 
side value was close to the one found by Charach [58]. 
Nevertheless, the melting temperature of PCM can be the 
arithmetic mean of HTF inlet temperature and temperature of 
surroundings for maximum exergy performance. As already stated, 
it is more appropriate to focus on complete cycle performance 
rather than performance during charging or discharging process 
alone. Hence, it can be mentioned that Eq. (15b) is more 
appropriate than Eq. (15a). 

The selection of appropriate PCM with the required optimum 
melting temperature is followed by selection of operating 
condition, i.e. the initial temperature for PCM. The initial 
temperature of PCM should be in accordance with the inlet 
temperature of HTF. As explained earlier, the inlet temperature of 
HTF should be as low as possible for charging and as high as 
possible for discharging. Hence, it is clear that the initial 
temperature for PCM should also be low for charging and high 
for discharging. This means the PCM should be at subcooled and 
superheated conditions, respectively, at the beginning of charging 
and discharging processes. Venkataramayya and Ramesh [27] have 
shown that the overall exergy efficiency of the system with 
subcooled and superheated PCM is higher than that of system with 
PCM kept at its melting point. However, the superiority of the 
former than the latter was found to be established only beyond the 
optimum melting point. This prompts to conclude that the initial 
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temperature of the PCM may be chosen in accordance with its 
melting temperature. 


5.3. Influence of system dimensions 


In any LHTS module, there are generally two compartments 
separated by a wall which is the heat transfer surface. One 
compartment is loaded with PCM and HTF flows through the other 
one. As an example, a typical shell and tube LHTS module showing 
PCM and HTF compartments is presented in Fig. 4. The volume of 
the PCM unit is determined by the quantity of PCM to be loaded, 
which in turn is determined by the quantity of heat to be stored or 
retrieved. However, for a fixed volume, the dimensions of the unit 
can still be adjusted to obtain maximum possible exergy 
performance. The dimensions of the PCM unit are expected to 
influence the heat transfer mechanism of the phase change 
process. Since the solidification of PCM is always governed by 
conduction heat transfer [60], it is generally not affected by the 
dimensions of the unit. On the other hand, the most part of the 
melting process is dominated by natural convection in the liquid 
PCM, although it is initially driven by conduction. The convection 
dominated melting process is reported by many investigators for 
almost all configurations. For example, in rectangular module by 
Stritih [61], in cylindrical module by Regin et al. [6] and in spherical 
module by Tan [62]. However, Shatikian et al. [63] have shown that 
the dominance of natural convection is significant only in larger 
rectangular systems. Similar results are presented by Ettouney 
et al. [64] for spherical modules. Hence, the increase in vertical 
dimension of the PCM unit, for instance, the diameter of the shell in 
case of shell and tube module, promotes the natural convection in 
the liquid PCM. The increase in natural convection may enhance 
the heat transfer rate. However, the higher rate of entropy 
generation is also expected due to faster interface motion and 
liquid viscosity as a result of enhanced natural convection. Erek 
and Dincer [33,51] have reported that increasing shell radius 
reduces the entropy generation number considerably. However, 
these works are limited to discharging process only. The effect of 
diameter/height of PCM unit on the overall exergy performance is 
not reported in the literature so far. 

It is understood that decreasing the length and increasing the 
diameter of the HTF unit could reduce the pressure drop in the HTF 
flow. This would reduce the entropy generation due to pressure 
drop irreversibilities. However, decrease in the length of the tube 
reduces the heat transfer surface area. According to El-Dessouky 
and Al-Juwayhel [52], entropy generation increases with decrease 
in heat transfer surface area. In fact, Charach and Zemel [55] have 
shown that number of transfer (NTU) which is a measure of heat 
transfer surface area (proportional to length) affects the entropy 
generation in two ways. They are 


e The entropy generation due to heat transfer irreversibilities 
decreases exponentially as the NTU increases. 


HTF out 


Fig. 4. PCM and HTF compartments in a typical shell and tube LHTS module. 
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Fig. 5. Arrangement of multiple PCMs in a shell and tube module. 


e The entropy generation due to pressure drop irreversibilities 
increases linearly as the NTU increases. 

Hence, there exists an optimum value for NTU, i.e. heat transfer area 

or length. This indicates that the length of the tube should be in such 

way that NTU of the unit is optimum. 


6. Exergy analysis with performance enhancement techniques 


In spite of various advantages like high energy density, nearly 
isothermal operation, etc., the phase change material loaded in the 
LHTS unit possesses a very low thermal conductivity, which 
drastically affects the performance of the unit. The effect of the 
lower value of conductivity is reflected during energy retrieval or 
withdrawal with an appreciable temperature drop during the 
process. As a result, the rate of phase change process (melting/ 
solidification of PCM) has not been up to the expected level and the 
large scale utilization of LHTS units, remains unsuccessful. 
Therefore, to tackle the above-mentioned drawbacks, it becomes 
necessary to improve the thermal performance of the LHTS units 
employing PCMs. Numerous investigations on various perform- 
ance enhancement techniques, both numerically and experimen- 
tally have been reported so far. The various techniques adopted for 
enhancing the thermal performance of the LHTS units are 
enumerated below: 


e Using extended surfaces [61,65-68]. 

e Employing multiple PCMs method [69-72]. 

e Thermal conductivity enhancement with additives. 
e Microencapsulation of PCMs [73-75]. 


For the recent review on the various performance enhancement 
techniques employed for LHTS units, readers are referred to the 
paper by Jegadheeswaran and Pohekar [76]. 

The studies on LHTS units employing performance enhance- 
ment techniques have mainly focused on evaluation of melting/ 
solidification time, heat transfer rate and amount of energy stored/ 
retrieved in comparison with those of system without enhance- 
ment techniques. From the results of these studies, it may be stated 
that all the proposed techniques have good potential in enhancing 
the thermal performance of LHTS systems. However, the enhance- 
ment due to employment of the above listed methods except 
multiple PCMs method is invariably assessed only through energy 
analysis. Even for multiple PCMs method, the exergy techniques 
are employed only by few investigators [38,40,47,77,78]. 


6.1. Multiple PCMs method 


As the name implies, the LHTS system is packed with more than 
one PCM of different melting temperatures for multiple PCMs 
method. In case of conventional single PCM system, the tempera- 
ture of HTF decreases in the direction of flow. As a consequence, the 
heat transfer rate decreases and thus, leads to poor performance of 


the unit. If multiple PCMs of different melting temperatures are 
packed in the unit in the decreasing order of their melting points, 
then nearly a constant temperature difference can be maintained 
all along the system, even though the HTF temperature decreases 
in the flow direction. This leads to almost a constant heat flux to the 
PCM. During discharging, if the HTF flow direction is reversed then 
the PCMs remain in the increasing order of their melting points. 
Once again nearly constant heat flux but from the PCM to HTF is 
possible. Recently, Bi et al. [79] have noticed that properly 
controlled uniform phase change rate would result in optimal 
operating characteristic of LHTS system, i.e. minimum entropy 
generation. Moreover, poor heat transfer rate from HTF to PCM 
during charging process results in relatively high temperature HTF 
at exit. Adebiyi [80] has shown that significant exergy loss would 
occur, if HTF exits the system at a rather high temperature for a 
long time during the storage process. This exergy loss can be 
expected to be less, if multiple PCMs are employed. This is because 
the HTF temperature at exit would be relatively less in multiple 
PCMs systems. The employment of multiple PCMs in shell and tube 
heat exchanger is illustrated in Fig. 5. 

The exergy analysis employed for LHTS systems using multiple 
PCMs may aim at the following: 


e To investigate the role of multiple PCMs in improving the exergy 
performance. 

e To identify the number of PCMs required. 

e To evaluate the optimum melting point distribution. 


As stated above, the use of multiple PCMs is basically for 
enhancing the charging or discharging rate. Watanabe and 
Kanzawa [38] have shown that rapid charging or discharging 
process as a result of multiple PCMs, led to high charging/ 
discharging/overall exergy efficiency under all operating and 
design conditions. Gong and Mujumdar [40] have found that a 
maximum of around 100% increase in overall exergy efficiency is 
possible with three PCMs. It is also reported that the increase can 
go up to 200% with five PCMs. Domanski and Fellah [77] 
investigated the improvement in the overall exergy efficiency 
when using two PCMs. The overall efficiency was calculated using 
entropy generation number and it was found that the increase was 
around 40% as compared to single PCM system. Similarly, the 
overall exergy efficiency of three PCMs system was found to be 74% 
higher than that of single PCM system by Gong and Mujumdar [78]. 
This was due to the fact that the charge/discharging process was 
around 37% faster in three PCMs system. Recently, Kousksou et al. 
[47] also noticed significant reduction in irreversibility when more 
number of PCMs were used. Hence, it is clear that multiple PCMs 
method is very promising from exergy point of view too. 

Since the exergy performance increases with increase in 
number of PCMs, it is also important to identify the number of 
PCMs that can be used. This important issue is addressed by Gong 
and Mujumdar [40]. The authors compared the overall exergy 
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efficiencies of the system using two, three and five PCMs with that 
of system using single PCM for different NTU values. The results 
revealed that significant improvement in exergy efficiency was 
obtained, when NTU were higher. On the other hand, the increase 
in exergy efficiency with multiple PCMs was marginal, when NTU 
values were smaller. This indicates that employment of more 
number of PCMs can be justified only if the heat transfer surface 
area is more. Hence, it can be stated that the number of PCMs needs 
to be chosen in accordance with the available heat transfer surface 
area. 

Since the PCMs are having different melting temperatures, 
attention should also be given to melting temperature distribu- 
tion. Domanski and Fellah [77] attempted numerically to obtain 
the melting points corresponding to the best exergy efficiency. 
However, this study is limited to two PCMs only. It is found that 
the optimum melting temperature of first PCM increases linearly 
with inlet temperature of the HTF. As it is known, this is to 
reduce the entropy generation due to the temperature difference 
between HTF and PCM. On the other hand, the second PCM (last 
one) should have a melting temperature close to temperature of 
atmosphere in order to minimize the amount of exergy loss to 
the surroundings. Watanabe and Kanzawa [38] have stated that 
the difference between the highest and lowest melting 
temperatures has enormous effects on the charging and 
discharging rates and thus, on the exergy performance. The 
authors have presented a simple expression for the optimum 
difference between the highest and lowest melting tempera- 
tures as follows. 

For charging, 


NTU 
T m( first PCM) oa T m(last PCM) a T+ (NTU/2) (THIF in Tmiave)) 


(16a) 


where Tm(ave) is the average melting temperature of the PCM 
located at the middle. 
For discharging, 


NTU 
Tiny first PCM) — Tm(last PCM) = T+ (NTU/2) (Tm(ave) — Ture,in) 


(16b) 


The two different conditions proposed for charging and 
discharging are difficult to satisfy as same set of PCMs are 
generally used for the cyclic operation. Moreover, when more than 
one PCM is involved, it is also necessary to examine the optimum 
conditions for melting temperatures of intermediate PCMs. This 
particular issue is not addressed by Watanabe and Kanzawa [38]. 
However, Gong and Mujumdar [40] obtained optimum melting 
temperatures for all PCMs, by maximizing the overall exergy 


Table 6 
Exergy based investigations on LHTS systems using multiple PCMs. 


efficiency of multiple PCMs module. The results of case studies 
based on this procedure indicate that the optimum melting 
temperatures of PCMs are approximately a geometrical regression 
i.e. Tm1/Tm2 & Tm2/Tm3 & Tm3/Tm4 S..  » -œ Tmn — 1/Tmn. Moreover, 
the ratio between the melting temperatures of two successive 
PCMs is found to be dependent on NTU of the system. 

It is clear from the above discussions that only Gong and 
Mujumdar [40] have investigated the optimum conditions for 
melting temperatures of all the PCMs in the system. The 
procedure adopted by the authors seems to be similar to that 
one used for single PCM system. However, the authors have 
reported that the maximization of overall efficiency of 
multiple PCMs system involves complex and tedious mathe- 
matical manipulation. The procedure becomes more and more 
complex as the number of PCMs increases. Hence, guidelines for 
optimum melting point distribution may not be drawn only 
based on this procedure. It is felt that comparative studies 
focusing on different proposals discussed above may lead to a 
concrete conclusion. It also requires more research as only a 
limited number of works is reported for exergy analysis of 
multiple PCMs system. Table 6 summarizes the outcomes of the 
exergy based investigations on LHTS systems using multiple 
PCMs. 


6.2. Addition of fins, thermal conductivity enhancement and 
microencapsulation of PCM 


Fins are employed for increasing the heat transfer area and thus, 
heat transfer rate. The conduction governed discharging (solidifi- 
cation) rate can be increased monotonically with increase in 
number of fins. On the other hand, there exists an optimum 
number of fins beyond which the convection dominated charging 
(melting) rate decreases [81]. This is because more number of fins 
dampens the natural convection in the liquid PCM. Hence, the 
addition of fins has direct impact on natural convection apart from 
increasing the heat transfer area (NTU). As we have seen, both 
natural convection in the liquid PCM and NTU can significantly 
affect the entropy generation. Therefore, the optimum number of 
fins should be computed based on exergy analysis. Since no work is 
reported on exergy analysis of finned system, it is not yet clear how 
the addition of fins enhances the exergy performance of LHTS 
systems. 

The conventional PCMs-both organic and inorganic possess 
very low thermal conductivity ranging from 0.1 to 0.6 W/m K. 
Hence, the additives of various forms as given below are used to 
improve the thermal conductivity. 


Ref. Number of PCMs Operation Approach (Exergy/ Operating/design Important findings 
mode Entropy generation) parameters considered 
[38] Seven Charging Entropy generation - Minimum irreversibility when melting 


HTF inlet temperature, NTU 


HTF inlet temperature, 


temperatures are linear 


Exergy efficiency increases by two to three times with 
multiple PCMs 

Optimum melting points of PCMs are approximately a 
geometric progression 


Faster charging or discharging leads to high exergy efficiency 


Initial temperature of PCM 


[40] Two, three and five Cycle Exergy 
[47] Fifteen or thirty Cycle Exergy 
[77] Two Cycle Exergy 
[78] Three Cycle Exergy 


Velocity and inlet 
temperature HTF, NTU 


Latent heat, NTU 


Melting point difference between first and last PCM is 
key for optimum exergy performance 


Melting temperatures of first and last PCMs to be close to 
HTF and atmospheric temperatures, respectively, for 
optimum exergy efficiency 

Increase in exergy efficiency due to multiple PCMs is not 
affected by Latent heat of PCM 
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e High conductivity porous material such as graphite [82-85]. 

e High conductivity and low density materials such as carbon 
fibers [86-88]. 

e High conductivity micro particles such as Aluminium, Copper, 
Silver, etc., [21,89-91]. 

e High conductivity metal structures [53,92,93]. 


In general, however, these additives may lead to loss of storage 
capacity of pure PCM as the mass of PCM decreases. Therefore, it is 
important to estimate the optimum mass/volume fraction of 
additives. It can be stated that the optimization of mass/volume 
fraction of additives should not be from storage point of view 
alone. This is because the additives of any form are expected to 
contribute to the total entropy generation. The same reason can be 
given for the need of exergy analysis in case of microencapsulation 
as solid structures are used for encapsulation of PCM. 


7. Thermoeconomic analysis 


It is important to note that the analyses of different researchers 
have firmly established the need to employ second law techniques 
to design thermodynamically efficient thermal energy storage 
systems. It is well known that, when a component of a thermal 
storage system displays a deterioration in terms of irreversibility, 
not only its performance, but also those of the remaining units 
which make up the system, are affected. In general, this makes 
difficult to understand what is the effect of a deterioration increase 
for the overall system, and, consequently, to determine for anyone 
of the components characterized by irreversibility, how the overall 
performance of the system would increase, if that irreversibility 
was eliminated. Thermoeconomics and its application to engi- 
neering energy systems can help to answer these questions. 

As stated above, thermoeconomic analysis combines economic 
and thermodynamic analysis by applying the concept of cost 
(originally an economic property) to exergy (an energetic property) 
(see Valero and Lozano [94]). Most analysts agree that exergy is the 
most adequate thermodynamic property to associate with cost since 
it contains information from the second law of thermodynamics and 
accounts for energy quality. Exergy based thermoeconomics 
methods are also referred to as “exergoeconomics”. 

Thermoeconomic analysis is based on the representation of a 
thermal system by using a thermoeconomic model. This is done by 
describing the system by means of a productive structure [95,96]. 
The productive structure (also “functional diagram” or “structural 
scheme”) is a graph which shows a set of relations defining the 
interaction among components themselves and the environment. 
Each subsystem has one or more entering flows that represent the 
resources or fuel and exiting flows (product). Due to irreversi- 
bilities, in every component, product is smaller than fuel. So, the 
amount of resources per unit of exergy increases throughout the 
system. This leads us to the concept of cost. 

Many thermal engineers have studied thermoeconomics or 
exergoeconomics, and various methodologies have been sug- 
gested. According to Tsatsaroni [97], thermoeconomics methods 
can be subdivided into two categories, those based on cost 
accounting and those based on optimization techniques. Cost 
accounting methods help to determine actual exergy cost and 
provide a rational basis for pricing, while optimization methods 
are used to find the optimum design or operating conditions [98]. 

The main feature of the thermoeconomics methods is that they 
propose a cost balance equation applying the exergetic unit cost to 
the exergy balance equation according to a specific principle [99]. 
As it was confirmed in related works [100,101], the cost structure 
of the overall system is mostly affected by the entropy production 
rate (irreversilbility) at each component and is dependent on the 
level of aggregation of the system. The level of aggregation 


provides a breakdown of the total irreversibility among the system 
components. The chosen level of aggregation will affect the 
conclusion of the analysis. 

In the scope of energy storage systems, to the best of the 
authors’ knowledge very few investigations using thermoeco- 
nomics methods have been reported [41,102]. Demirel and Ozturk 
[41] have performed the cost accounting method to study a 
seasonal latent storage system for heating a greenhouse. In order to 
build up the cost formation in the system, structural diagrams have 
been used. To simplify the cost balance model, all the work 
interactions for fans, valves and pressure changes in each 
component are assumed to be negligible. The cost rate balance 
is established for both charging and discharging modes. They 
found that the feasibility of overall design and operation of the 
thermal storage system can be enhanced using the cashflow- 
diagram based on exergy considerations. Badar et al. [103] and 
Badar and Zubair [104] discussed the second law-based thermo- 
economic analysis of a sensible heat storage system. They 
described the procedure for optimizing the performance of a 
thermal energy storage system, which has no removal cycle. They 
also compared the results with that obtained from Bejan’s analysis 
[23]. Recently, Zubair and Al-Naglah [105] extended thermo- 
economic analysis of a storage system that has both energy storage 
and removal processes with convection by ideal gas. They found 
that the unit cost of irreversible losses due to heat transfer is 
relatively more sensitive to the number of units transfer when 
compared to that of due to pressure drop. To optimize the 
performance of the studied system they have introduced a new 
performance criterion described as the cost rate number. 

In general, the objective of thermoeconomic optimization is to 
minimize the annualized total cost of owning and operating the 
system. The total cost is composed of additive terms that measure 
the costs associated with the lost exergy and the annualized capital 
costs of the equipment. However, thermal energy storage (TES)- 
based systems are usually economically justifiable when the 
annualized capital and operating costs are less than those costs for 
primary generating equipment supplying the same service loads 
and periods. TES is mainly installed to lower the initial costs of the 
other plant components and operating costs. Lower initial 
equipment costs are usually obtained when large durations occur 
between periods of demand. Operating cost savings and the net 
overall costs should be assessed using life cycle costing or other 
suitable methods to determine which system is the most 
beneficial. From this analysis, the aim of thermoeconomic 
optimization based on exergetic analysis must be directed towards 
minimizing the exergetic cost of the system. So, the problem of 
thermoeconomics analysis applied to energy storage systems can 
be formulated as follows: Given a system whose limits have been 
defined and a level of aggregation that specifies the components 
which constitute it, how to obtain the cost of all the flows that 
become interrelated in this system. This statement is a cornerstone 
in thermoeconomics. It is also necessary to mention that the 
results of the thermoeconomic analysis of the energy storage 
systems will mainly dependent upon: 


- the storage period required (i.e., diurnal, weekly or seasonal), 

- the defined limits of the system and the definition of resources 
are always partial and refer to the system under study. Generally 
speaking, energy, raw materials, PCMs, economy and labour 
resources are those which are put at the disposal of the system 
within its limits of analysis and with known unit prices. 


8. Conclusions 


In the present paper, the various techniques adopted for the 
exergy based performance evaluation of LHTS units and subse- 
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quent optimization of the system have been reviewed. This review 
intends to motivate the researchers to focus more on exergy 
analysis in the future research works on LHTS systems. To 
conclude, the following points are highlighted. 


e Exergy based performance evaluation of LHTS system can be 
stated as more perspective measure than energy based one as it 
reflects the true potential of the system and the economic value 
of the storage/recovery operation. 

Exergy analysis of any form should consider the complete cycle 
operation of LHTS systems rather than considering either storage 
or recovery mode alone. 

Optimization of any LHTS module requires rigorous exergy 
analysis along with energy analysis. 

Entropy generation analysis has a significant importance 
especially when it comes to optimization of design and operating 
parameters of LHTS systems. 

Studies still remain on large scale especially for the comparative 
evaluation of various performance enhancement techniques 
employed for LHTS systems. 

The installation and operation of any LHTS system can be 
justified only if the exergetic cost of the system is minimum and 
exergy based thermoeconomic optimization is of greater help in 
minimizing the exergetic cost of the system. 
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